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ABSTRACT 

We isolate samples of 733 bright {g < 18) and 437 faint (g > 18) high-Galactic 
latitude blue horizontal branch stars with photometry and spectroscopy in the Sloan 
Digital Sky Survey (SDSS). Comparison of independent photometric and spectroscopic 
selection criteria indicates that contamination from F and blue-straggler stars is less 
than 10% for bright stars (g < 18) and about 25% for faint stars (g > 18), and this 
is qualitatively confirmed by proper motions based on the USNO-A catalog as first 
epoch. Analysis of repeated observations shows that the errors in radial velocity are 
~ 26kms~ 1 . 

A relation between absolute magnitude and color is established using the horizontal 
branches of halo globular clusters observed by SDSS. Bolometric corrections and colors 
are synthesized in the SDSS filters from model spectra. The redder stars agree well in 
absolute magitude with accepted values for RR Lyrae stars. The resulting photometric 
distances are accurate to about 0.2 magnitudes, with a median of about 25 kpc. Modest 
clumps in phase space exist and are consistent with the previously reported tidal stream 
of the Sagittarius dwarf galaxy. 

The sample is tabulated in electronic form in the online version of this article, or by 
request to the authors. 

Subject headings: Galaxy: halo — Galaxy: structure — stars: horizontal branch 

1. Introduction 

Blue horizontal branch (BHB) stars have long been important objects of study as tracers of the 
Galactic potential (Pier 1983; Sommer-Larsen, Christensen, & Carter 1989; Clewley et al. 2002). 
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Their constant absolute magnitude allows accurate space positions to be determined. Historically, 
many researchers have used samples of BHB stars, or other standard-candle stars such as RR 
Lyraes, to study the local halo, i.e. the vicinity of the solar neighborhood where typical distances 
are of order the solar Galactocentric radius Ro = 8 kpc. These samples can be used to constrain the 
kinematics of the local halo (Pier 1984; Layden et al. 1996). The distant halo (i.e. typical distances 
3> Ro), however, is less well constrained (but see Sommer-Larsen et al. (1997)). Isolating a sample 
of distant BHB stars would enable study of the structure and kinematics of the Galactic halo on 
larger scales (Sirko et al. 2003, hereafter Paper II), and estimates of the mass of the Galaxy. 

The Sloan Digital Sky Survey (SDSS), currently halfway through a photometric and spectro- 
scopic survey of a quarter of the sky, offers an excellent opportunity to isolate and study BHBs. 
The SDSS spectroscopic survey targets objects between about g = 15.5 to 19 (the bright limit is set 
to avoid contamination in adjacent spectroscopic fibers at the slit end of the fiber bundle), although 
a significant number of stars outside this range is also observed. Thus BHB stars at distances from 
~ 7 to ~ 60 kpc from the sun appear in the spectroscopic survey. 

This paper presents 1170 BHB stars selected from SDSS that probe the outer halo. The plan 
of the paper is as follows. Pertinent characteristics of the Sloan Digital Sky Survey are summarized 
in §2. §3 describes the preliminary color cut. §4 discusses spectroscopic BHB classification criteria. 
For the faint stars, the spectroscopic classification breaks down so we resort to a more stringent 
color cut described in §5. Radial velocity errors are evaluated in §6. §7 discusses the calculation of 
bolometric corrections for accurate distances, and evaluates magnitude errors. Some sanity checks 
are presented in §8, including an analysis of the limited proper motion data available, and a map 
of the Galactic halo which shows the conspicuous Sagittarius tidal stream (figure 13). The results 
are summarized in §9. The sample is available in Table 3 in the electronic version of this paper. 

2. The Sloan Digital Sky Survey 

The Sloan Digital Sky Survey (SDSS: York et al. 2000) is a project to image about 1/4 of 
the sky at high Galactic latitude to a depth of about 21 to 22 magnitudes in five broad bands 
covering the entire optical range and to obtain spectra of about a million galaxies and 100,000 
quasars selected from the imaging data. While galaxies and quasars are the primary targets of 
the spectroscopic survey, the SDSS also obtains a vast amount of stellar data. Stellar spectra are 
acquired by several routes: observations of calibration stars, observations of various classes of star 
to be used as "filler" observations when a given region of sky has too low a density of primary 
targets, and via objects targeted as candidate galaxies and quasars that turn out to be stars. 

The SDSS uses a dedicated 2.5 meter telescope and a large format CCD camera (Gunn, Carr, 
Rockosi, & Sekiguchi et al. 1998) at the Apache Point Observatory in New Mexico to obtain images 
almost simultaneously in five broad bands (u, g, r, i, z) centered at 3551, 4686, 6166, 7480, and 8932 
A respectively (Fukugita et al. 1996; Gunn, Carr, Rockosi, & Sekiguchi et al. 1998). The imaging 
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data are automatically processed through a series of software pipelines which find and measure 
objects and provide photometric and astrometric calibrations to produce a catalogue of objects 
with calibrated magnitudes, positions and structure information. The photometric pipeline Photo 
(Lupton et al. 2001, and in preparation) detects the objects, matches the data from the five filters, 
and measures instrumental fluxes, positions and shape parameters. The last allow the classification 
of objects as "point source" (compatible with the point spread function [psf]), or "extended." The 
instrumental fluxes are calibrated via a network of primary and secondary stellar flux standards 
to AB,, magnitudes (Fukugita et al. 1996; Hogg et al. 2001; Smith et al. 2002). We use the "psf 
magnitude," which represents the magnitude of the best-fit psf to the object, as determined by 
Photo version 5.4.25 (2003 June). The psf magnitudes are currently accurate to about 2% in g,r, 
and i and 3-5% in u and z for bright (< 20 mag) point sources. The SDSS is 50% complete for point 
sources to (u, g, r, i, z) = (22.5, 23.2, 22.6, 21.9, 20.8), and the full- width at half-maximum (FWHM) 
of the psf is about 1.5" (Abazajian et al. 2003). The data are saturated at about 14 mag in g, r 
and i and about 12 mag in u and z. Astrometric calibration is carried out using a second set of 
less sensitive CCDs in the camera, which allows the transfer of astrometric catalogue positions to 
the fainter objects detected by the main camera. Absolute positions are accurate to better than 
0.1" in each coordinate (Pier et al. 2003). 

The SDSS spectroscopic survey is carried out by targeting objects selected from the imaging 
catalogues according to pre-determined and carefully tested criteria (Stoughton et al. 2002). The 
targeted stars, quasars and galaxies have limiting magnitudes of around 18 — 20 mag (Richards 
et al. 2002; Eisenstein et al. 2001; Strauss et al. 2002). BHB stars appear in the SDSS spectroscopic 
sample via two main routes. At redshifts around 2.7, the quasar color locus crosses the BHB color 
locus (Richards et al. 2002) and many BHB stars are targeted as candidate quasars. BHB stars are 
also photometrically selected to serve as filler observations in regions of the sky with a low density 
of primary targets. The sky locations of the target objects are mapped to a series of plug plates 
which feed a pair of double spectrographs (A. Uomoto, S. Smee et al. unpublished) via optical 
fibers. The spectrographs jointly observe 640 objects (including sky positions and standard stars) 
simultaneously, and provide continuous wavelength coverage from about 3800 A to 9500 A at a 
resolution of about 1800 — 2100. The nominal total exposure time for each spectroscopic plate is 45 
minutes (15 minutes x 3 exposures), plus observations of flat field and arc standard lamps for every 
plate. The entrance aperture of the fibers is 3". The spectra are matched to the photometric objects 
using a video mapping system for the plug plates (D. Schlegel and D. Finkbeiner, unpublished). 
The spectra are recorded on CCD detectors. The data are optimally extracted, the sky spectrum is 
subtracted, and wavelength, flat-field and spectrophotometric calibrations are applied automatically 
by the Spectro-2D software pipeline, written by D. Schlegel and S. Buries (unpublished). 

Each spectroscopic plate is a circle of diameter 3° on the sky. The present work is based on the 
data of 937 unique spectroscopic plates taken to date, through 2003 Aug 6. In addition, 77 plates 
were observed redundantly, but in defining the BHB sample we exclude redundant observations by 
adopting the plate with the highest S/N. However, we use the information from these redundant 
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observations in §6 below to evaluate velocity errors. 

The radial velocity is found by fitting the spectrum to template spectra, using an automated 
fitting program (Spectro_Brightstars) written by D. Schlegel (unpublished). The basic procedure 
is as follows. Prugniel & Soubiran (2001) have published a catalogue of about a thousand spectra 
of stars with a wide range of properties (spectral types, metallicities etc.) observed at high spectral 
resolution with the ELODIE spectrograph. One version of the ELODIE spectra is binned to R 
~ 20,000 and spectrophotometrically calibrated to about 5%. Brightstars was used to select 
template spectra from the SDSS data which match the ELODIE spectra. These template spectra 
are typically composites from many SDSS spectra, and are set up to provide relatively coarse 
sampling of the spectral type grid. The template spectra are aligned to zero heliocentric velocity. 
Each SDSS spectrum is then compared with the template spectra in wavelength space, and the 
template and wavelength shift that give the best match is automatically measured. The stars in 
the present sample were almost all best-fit with spectral types between B6 and AO. The radial 
velocities are corrected to the heliocentric standard of rest, and throughout this paper are denoted 



3. Color cut 

BHB stars have very characteristic colors in the high-Galactic latitude stellar sample observed 
by the SDSS imaging survey, and photometric selection of candidate objects is straightforward. 
BHB stars are bluer in g — r than most halo stars because the halo main sequence turnoff lies at 
spectral type F/G. The large Balmer jump in BHB low-gravity stars gives them redder u — g colors 
than other blue objects such as low redshift quasars and white dwarfs, and indeed the SDSS u and 
g filters are specifically designed to provide optimal photometric separation of low-redshift quasars 
(Gunn, Carr, Rockosi, & Sekiguchi et al. 1998; York et al. 2000; Richards et al. 2002). The SDSS 
g — r vs. u — g color-color diagram in figure 1 shows all objects in the spectroscopic sample of 937 
plates used herein that were subsequently identified as stars by Brightstars, which is accurate 
enough that quasars and galaxies are virtually nonexistent in the sample. (Throughout this paper 
all magnitudes are corrected for Galactic extinction using the maps of Schlegel, Finkbeiner, &: Davis 
(1998).) The region occupied by BHB stars is indicated with a "color cut box," with boundaries 
0.8 < u — g < 1.35, —0.4 < g — r < 0.0. The boundaries of this selection box are similar to those 
used by Yanny & Newberg et al. (2000), and figure 8 (see §5) shows that our final sample of bright 
BHB stars is well contained within the boundaries. Stars that lie in this region of color-color space 
pass what we will refer to as the "color cut." 

Although BHB stars are effectively selected through the color cut described above, there re- 
mains a significant degree of contamination from two other types of star. The locus of cooler F stars 
is actually quite outside of the color cut box, to the red in g — r (Yanny & Newberg et al. 2000), but 
there are so many F stars that they leak into the BHB candidate sample, via intrinsic variations 
and photometric errors, with numbers comparable to the number of BHB stars. To identify these 
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cooler stars merely requires a measure of temperature independent of u — g or g — r. Somewhat 
more challenging to identify are the A-type stars with higher surface gravity (Preston & Sneden 
2000; Yanny & Newberg et al. 2000; Clewley et al. 2002). These are probably field blue stragglers, 
which are some 2 magnitudes less luminous (Preston & Sneden 2000). While the Balmer jump as 
measured by u — g can separate out many of these stars, some remain in the sample. They can be 
removed using their broader Balmer lines as described in the next section. 

4. Balmer line analysis cuts 

The upper panel of figure 2 presents a typical BHB spectrum, showing the deep Balmer lines 
of an A-type star. In this section we describe two methods for analyzing these lines: the widely 
used .Do. 2 method (Pier 1983; Sommer-Larsen & Christensen 1986; Arnold & Gilmore 1992; Flynn, 
Sommer-Larsen, &i Christensen 1994; Kinman, Suntzeff, & Kraft 1994; Wilhelm, Beers, & Gray 

1999) and the scalewidth-shape method (Clewley et al. 2002). We focus on the H7 and H<5 lines 
because the the H/3 and Ha lines have lower continuum levels, and the lines above H<5 are much 
more closely spaced, making the determination of the continuum difficult (Yanny & Newberg et al. 

2000) . 

To analyze the lines, the continuum must first be fitted and divided out. For the -Do. 2 method, 
following Yanny & Newberg et al. (2000), we take the continuum to range from 4000 to 4500 A 
with 60 A masks covering the positions of the H7 and H<5 lines. We then divide the spectrum by 
the best-fit sixth-order Legendre polynomial to this extracted continuum. Finally, the spectrum 
is boxcar-smoothed with a boxcar length of 5 pixels. An example normalized and extracted SED 
is shown in the lower panel of figure 2. Similarly, for the scalewidth-shape method, we take the 
continuum extraction rule of Clewley et al. (2002) (which ranges from 3863 to 4494 A) and fit to a 
fourth-order Legendre polynomial. After the continuum is divided out with this prescription, the 
SED looks very similar to the example SED in the lower panel of figure 2. 

The following terminology is used in this section and throughout. A "method" refers to a 
general algorithm, which may or may not be desirable for isolated use, which employs one or more 
"criteria" for discriminating BHB stars. The -Do. 2 method and the scalewidth-shape method are 
presented in the following. A "cut" is a method or combination of methods that we actually use 
to eliminate contaminants from the sample. The combination cut discussed below incorporates the 
-D0.2 and scalewidth-shape methods. 

4.1. The D , 2 method 

Since BHB stars have low surface gravity, their Balmer lines are narrower than the smaller main 
sequence blue stragglers. The -Do. 2 method discriminates BHB stars from BS stars by determining 
the value of -Do. 2, the width in angstroms of the Balmer line at 80% of the continuum (Yanny & 
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Newberg et al. 2000). Figure 3 demonstrates how this parameter can discriminate a BHB from a 
BS star. Although simple, this method has been found to be fairly robust and can provide a better 
measurement of the linewidth than other more complicated procedures (e.g. Flynn, Sommer-Larsen, 
& Christensen 1994). 

The -Do. 2 method nominally only applies to objects with A-type spectra. However, as men- 
tioned in §3, a large contaminant of our sample is F stars. Since F-type stars have weaker Balmer 
lines than A-type stars, one way of measuring temperature is by measuring the depth of the H7 or 
H5 line. Here we define / m as the flux relative to the continuum at the minimum of the line. Thus, 
higher values of / m for the Balmer lines signify cooler stars. A plot of Z?o.2 vs. fm f° r stars passing 
the color cut is shown in figure 4. In this figure, only stars with g magnitudes < 18 are shown, so 
that the trends in the plot are not obscured by the noise in the data of fainter stars (see figure 5, 
to be discussed in §4.3). 

The concentration of stars near (/mj-Do.2) = (-30, 26 A) represents the BHB stars, and those 
stars with larger values of Do. 2 are BS stars. Figure 4 also suggests that the locii of BS and BHB 
stars propagate to larger / m and smaller £>o.2 until the two locii converge at about (/mj-Do.2) = 
(.4, 23 A). In other words, for a given .3 < / m < .4, the distribution of stars in -Do. 2 is bimodal, 
but the stars are most cleanly separated for the smaller values of f m and not separated at all 
when / m ~ .4. Even though figure 4 doesn't show the trend clearly, we return to this point in 
§4.3 and show that the stars with smaller values of -Do. 2 but with .34 < / m < .4 are the other 
horizontal branch stars; RR Lyraes, for example, fall into this region of the horizontal branch. 
Because horizontal branch stars are evidently more difficult to discriminate from main sequence 
stars at lower temperatures, we adopt the criteria / m < .34, -Do. 2 < 28 for likely BHB stars for the 
He) -D0.2 method. 



4.2. The scalewidth-shape method 

Clewley et al. (2002) have proposed the scalewidth-shape method for discriminating BHB stars 
based on the Sersic profile for Balmer lines (Sersic 1968): 



y = 1.0 — aexp 



(1) 



Here, y is the normalized flux and x is the wavelength in A. We fit the normalized extracted 
spectrum (where the bounds of extraction are given by Clewley et al. (2002)) to the Sersic profile 
with three free parameters: a, b and c. The parameter xq is assumed to be the nominal location 
of the Balmer line corrected for redshift, which has already been determined. After fitting both 
Balmer lines of every star to the three-parameter (a, b, c) Sersic profile, there is a reasonably small 
dispersion in the determination of a: for the H*y line the average is .690 and the root-mean-square 
(rms) is .019; for H5 the average is .740 and the rms is .018. We fix a at its average value among 
stars passing the BHB criteria of the H5 -Do. 2 method, which is not objectionable because the .Do. 2 
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method is independent of the scalewidth-shape method. The Balmer lines are then refitted to the 
two- parameter (6, c) Sersic profile. This refitting is important because errors in a are correlated 
with the errors in b and c (Clewley et al. 2002). For this work we adopt a 7 = .690 for the H7 line 
and as = .740 for the H<5 line. The optimum value of a most likely depends on the resolution of 
spectroscopy and similar factors. 

Fitting to the Sersic profile thus provides the two parameters b and c for both the H7 and H5 
lines, for every star. Clewley et al. (2002) show that BHB and BS stars separate rather cleanly in 
a plot of b against c. Furthermore, they show that c is a measure of temperature: cooler stars have 
smaller values of c. 

4.3. Determination of criteria for the Do. 2 and scalewidth-shape methods 

We concentrate on the set of brighter stars with g < 18 because the noise in the spectra 
of fainter stars throws large errors in the parameter determinations for the .Do. 2 method and the 
scalewidth-shape method. This point is illustrated in figure 5. We will return to the selection of 
BHBs from the fainter star sample in §5. 

Figure 6 shows the parameters of the sample of BHB candidates with g < 18 as determined by 
the two linewidth analysis methods described in the preceding sections. Each star is represented 
by the same color in all four panels. In this figure the H<5 Do. 2 method is represented in the lower 
left panel, and stars which pass the BHB criteria for this method are colored blue, stars which are 
likely blue stragglers are green, and all others (cooler stars) are red. It is readily apparent that the 
four independent parameter spaces map BHB stars, BS stars, and cool stars to different regions. 

The exact BHB criteria for each method were determined through a combination of adjusting 
the bounds of the criteria by eye and comparing the results of one method to the others. For 
instance, the scalewidth-shape method separates BHB and BS stars cleanly enough that there is a 
natural gap between the two distributions in a plot of b vs. c, as shown in figure 6. A parabola is 
qualitatively drawn to divide the two distributions. Similarly, with regard to the Do. 2 method, the 
distribution in Do. 2 of the stars with low values of / m naturally separates into the two populations 
of BHB and BS stars. The issue of where to draw the temperature criterion (/ m for the -Do. 2 
method and c for the scalewidth-shape method) is more sensitive because the "natural" gap is 
harder to see in the .Do. 2 method and seems to be nonexistent in the scalewidth-shape method. 
In fact, the upper right panel of figure 6 hints that the division between stars with low surface 
gravity and high surface gravity continues even for much cooler stars than BHB stars. The cooler 
population (0.8 < c 7 < 1.0) with lower surface gravity (lower & 7 ) is thus the horizontal branch 
stars which are not blue, e.g. RR Lyraes. This b vs. c plot corroborates the evidence in §4.1 that 
non-blue horizontal branch stars form a locus in the (/ m) -Do.2) plane which extends out to cooler 
temperatures corresponding to f m ~ .4. Figure 7 shows that the stars of the population that has 
smaller values of 6 7 also have smaller values of Do. 2- With this information, we deduce that the BHB 
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stars are located in the clump in the H<5 / m — D0.2 plane centered on (/ m , Do.2)h5 ~ (-30, 26 A) and 
we set the temperature criterion at / m = .34. A subsequent inspection of figure 6 reveals that most 
BHB stars have c 7 > f .0, so we adopt this as the temperature criterion for the H7 scalewidth-shape 
method. Table 1 lists the adopted BHB selection criteria for each method. 

4.4. The combination cut 

We define the combination cut as those stars which pass the criteria for (the H7 Do. 2 and 
H5 scalewidth-shape methods) OR (the H5 Do. 2 and H7 scalewidth-shape methods). We assume 
that the H7 Do. 2 and H<5 scalewidth-shape methods are completely independent of each other, as 
are the H<5 Do. 2 and H7 scalewidth-shape methods. This assumption is reasonable because no 
information is shared in the analyses of either line, and the Do. 2 and scalewidth-shape methods 
measure fundamentally different properties of the spectral energy distribution. Under this assump- 
tion, the contamination fraction for the combination cut would be ft C ombo ~ 2k 2 , where k is the 
contamination fraction of a single method. In Appendix A we show that the contamination fraction 
of any one of the four methods is probably k < 10%, so K C ombo ^ 2% under these assumptions. 
However, to be more conservative, the fact that k ~ 10% for all four individual methods suggests 
that K com bo is probably 10% or less; in Paper II we use K C ombo = 10% to study errors introduced 
by contamination into the kinematic analysis . 

5. The fainter stars 

As shown in figure 5, the linewidth analysis methods cannot be expected to perform accurately 
for fainter stars than g > 18; in fact the number of g > 18 stars that are identified differently in any 
two methods (~ 100 — 400) is of the same order as the number of g > 18 stars identified as BHB 
stars in one method (~ 200 — 500). Because these fainter stars have noisier spectra, not only are 
the algorithmic BHB criteria unreliable, but eye-inspection would be difficult at best. However, an 
inspection of the results of the combination cut applied to the bright {g < 18) stars, as shown in 
figure 8, reveals that BHB stars, BS stars, and cooler stars occupy different regions of color-color 
space, although there is significant overlap. In this figure BHB stars are defined as those which 
pass the combination cut; (green) BS stars are defined as those which do not pass the combination 
cut but lie in the BS region for the H<5 Do. 2 method; and red stars are everything else. The thick 
curved line with periodic circles will be explained in §7. 

The bounds of this particular plot are the same as the color cut discussed in §3, so it is evident 
that the BHB stars are more likely to be found in the region enclosed by the "piano" shape indicated 
than outside of it. Since SDSS photometry can be expected to be accurate to a fainter magnitude 
(g ~ 20, York et al. 2000) than can the combination method (see figure 5), we select BHB stars 
with g > 18 solely on the basis of color. The new color cut is indicated by the piano shape in 
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figure 8 (and is given by the region 0.85 < u — g < 1.3, —0.31 < g — r < —.13 excluding the elliptical 
subregion ((it - g - 0.85)/0.31) 2 + {{g - r + 0.13)/0.11) 2 < 1) and will be referred to here as the 
"stringent" color cut. 

What are the contamination and incompleteness for the stringent color cut with no spectro- 
scopic analysis? Assuming that the combination cut is 100% accurate for g < 18, there are 733 
BHB stars, 172 of which do not pass the stringent color cut. There are 199 stars that pass the 
stringent color cut that are not classified as BHB stars by the combination method. Most of these 
(133) are classified as blue stragglers by the H<5 Do. 2 method. Thus the contamination for the 
stringent color cut alone is 199/(733 - 172 + 199) » 25%. 

6. Duplicate stars and radial velocity errors 

Since the spectroscopic plates overlap on the sky, sometimes objects are targeted for spec- 
troscopy on two different plates. Of the 4515 stars passing the color cut (from unique plates), 81 
are duplicate observations of the same star on different plates. We remove these 81 redundant 
observations from our BHB sample. In addition, from the data of the 77 non-unique plates (see 
§2), there are 410 stars passing the color cut which duplicate observations of stars from the unique 
plates. Each spectroscopic observation of a star provides an independent measurement of radial 
velocity v r and, for the bright (g < 18) stars, an independent classification as a BHB star or non- 
BHB star. Figure 9 shows the 491 stars which have multiple spectra and therefore independent 
determinations of v r . The ordinate axis shows the difference in the radial velocity determinations 
Av r between each duplicate star and its corresponding star in our purified sample of unique stars. 
Stars which are classified as BHBs for both observations are colored blue, stars which are classified 
as non-BHBs for both observations are colored red, and stars which have conflicting classifications 
between the two observations are colored green. Brighter than g = 18, there are 91 stars which are 
doubly classified as BHBs, 131 stars which are doubly classified as non-BHBs, and 13 stars which 
have conflicting classifications. This is consistent with the 10% contamination determined in Ap- 
pendix A, assuming that there is about an equal number of contaminants in the BHB sample as the 
number of BHB stars which are not classified as such (i.e., assuming k ~ rj). Figure 9 also shows the 
dependence of velocity errors on magnitude; as expected, fainter stars tend to have greater velocity 
errors, presumably due to noiser spectra. The rank correlation coefficient of (g, |Aiy|) is 0.263. The 
velocity differences determined this way turn out to be consistent with the errors derived from the 
spectral fits for only about half of the stars, independent of the magnitude or whether or not a star 
is a BHB or BS star. This comparison implies that the errors have been underestimated and this 
is currently under investigation. The rms of Av r in figure 9 is 36.5 km s . Clearly the errors in 
radial velocity measurements are not gaussian, but a meaningful approximate result may still be 
obtained if they are assumed to be so. Since the distribution of Av r is based on the errors of two 
measurements, the error of one measurement is approximately Av r ~ 36.5/y2 ~ 26kms~ 1 . 
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7. The absolute magnitude of BHB stars 

In order to refine our photometric distance estimates, we have computed theoretical colors and 
bolometric corrections for BHBs. We integrate the latest Kurucz model atmospheres (downloaded 
from www . kurucz . harvard . edu) against the SDSS filter response functions at a standard airmass 
of 1.2 atmospheres, in the manner described by Lenz et al. (1998). Unfortunately, the latter authors 
report only colors, so it was necessary to repeat their calculations in order to obtain the bolometric 
corrections. Apart from wavelength, the Kurucz spectra are functions of three parameters: effective 
temperature T e g, surface gravity g = GM/r 2 , and metallicity relative to solar [M/H]. 

To derive absolute magnitudes and bolometric corrections, it is necessary to assume a relation- 
ship between the luminosity L and T e g. The models of Dorman, Rood, & O'Connell (1993) for the 
theoretical zero-age horizontal branch (ZAHB) at [M/H] < —1.5 can be summarized roughly by 

log(L/L ) « 1.58 - 0.73 log(T cfr /10 4 K) (2) 

over the range 8000 K < T c g- < 12000 K. From a detailed study of the horizontal branch of M5, 
however, Baev, Markov, & Spassova (2001) find that luminosity varies more weakly with T e g; in 
fact their results appear to be consistent with constant luminosity and mass, 

M&0.6M e , log— « 1.65 (3) 

Lq 

with a scatter of no more than ±O.O5M and ±0.05 dex, respectively, over the same temperature 
range. We adopt the mean values in equation (3), so that log g « 3.52 + 41og(T e fj/10 4 K). 

To give an idea of the theoretical uncertainties, we note that if in fact the relation (2) holds 
for our HB stars, then by adopting (3) we have overestimated their distance moduli by an average 
of 0.18 mag over the indicated range of T e g , and by neglecting the variation of L with T e g, we have 
increased the scatter in distance modulus by ~ 0.09 mag, which should be taken in quadrature with 
the photometric measurement error. 

The results of this exercise are shown in Table 2. Evidently, for the colors of interest to us, 
T cff G (8000, 12000) K, and the influence of metallicity is slight if [M/H] < -1. 

The sequence of theoretical horizontal branch stars, parametrized by temperature, is projected 
onto the (u — g, g — r) plane in figure 8. The five black circles on the thick line denote, starting at the 
top and progressing down and to the left, stars with absolute magnitudes g = (0.6, 0.55, 0.6, 0.7, 0.8). 
For every observed BHB star in our sample (colored blue in the figure) we define the most probable 
absolute magnitude associated with its (u— g, g — r) colors by simply finding the absolute magnitude 
of the point on the theoretical track that is closest to the observed star in this color-color space. 

The stringent color cut (piano shape in figure 8, see §5) was defined before the theoretical track 
was calculated, so the fact that the theoretical track satisfies the stringent color cut is a meaningful 
consistency check. 
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As a further consistency check, we have isolated several Galactic globular clusters observed 
by the SDSS photometric survey to date and plotted their color-magnitude diagrams (CMDs). 
Pal 3, NGC 5904, and NGC 6205 do not show a discernable horizontal branch, so they are not 
considered here. The g vs. g — r CMDs of Pal 5 (Yanny & Newberg et al. 2000; Ivezic et al. 
2000) and NGC 2419 are presented in figure 10, where the distance modulus, derived from cluster 
distances of d = 84.2 and 23.2 kpc for NGC 2419 and Pal 5, respectively (Harris 1996), has been 
subtracted from the dereddened apparent magnitude g. Stars within 8' of the center of NGC 2419 
and within 14' of the center of Pal 5 were selected. The SDSS automated photometry avoids 
analyzing very crowded fields, so the stars presented here, though within the stated radii, tend 
to lie in the outskirts of the clusters. Vertical lines at g — r = (—0.4, 0.0) indicate our color cut. 
The theoretical horizontal branch sequence is also plotted in both panels. The horizontal branch 
stars of NGC 2419 appear to be systematically fainter than the theoretical prediction, while those 
of Pal 5 seem to be systematically brighter; however, these systematic differences could easily be 
due to inaccurate distances d. By selecting from these two globular clusters the stars passing the 
color cut (0.8 < u — g < 1.35, —0.4 < g — r < 0.0) and lying in the region of the CMD for the 
horizontal branch (0.0 < g — 5 log (d/10 pc) < 1.5), we compare the directly observed absolute 
magnitude g — 5 log (d/10 pc) to the indirectly obtained theoretical absolute magnitude based on 
the (u — g,g — r) colors. For NGC 2419, the rms of the differences of the two among the 94 stars 
is .134 magnitudes, and the root-sum-square (rss) is .167 magnitudes 6 . For Pal 5, the rms is .126 
magnitudes and the rss is .202 magnitudes. The difference between the rms and rss is due to the 
nonzero mean, or systematic error as discussed above. Therefore, the magnitude error of a star, 
with absolute magnitude determined by comparison to the theoretical horizontal branch sequence 
described above, is of order .2 magnitudes. The error may be as small as ~ .13 magnitudes if 
the only cause of the systematic difference between directly observed absolute magnitudes and our 
indirectly obtained absolute magnitudes is incorrect distances for Pal 5 and NGC 2419. 



8. Sanity checks 

As figure 8 shows, BHB stars determined by the combination cut are in general redder in u — g 
and bluer in g — r than other objects such as blue stragglers. This result is consistent with the 
work of Yanny & Newberg et al. (2000), so it is reassuring. In this section we present several more 
sanity checks which verify that our sample of BHB stars is probably accurate. 



rms = ^^E(A< 7l -(Ap)) 2 , (A 5 > = rss = ±XAg? 
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8.1. Metallicity from the Call K line 

BHB stars are halo stars, and as such they should have lower metallicity than the disk star 
contaminants. The most promising spectral feature to use for a metallicity determination is the Call 
K line at 3933.7 A (the H line, at 3968.5 A, is not used because it is too close to He). Unfortunately, 
it is impractical to use the Call K line as a reliable indicator of metallicity in individual SDSS 
spectra because much higher spectral resolution would be necessary. However, the following exercise 
is executed anyway as a sanity check. By employing a continuum-subtraction algorithm very 
similar to that described in §4.1 for the H7 and H<5 Balmer lines, we find the value of / m for the 
Call K line; smaller values indicate larger metallicity. Figure 11 presents the distributions of this 
parameter for our sample of BHB stars and non-BHB stars passing the color cut. The figure shows 
different distributions for bright (g < 18) and faint {g > 18) stars because of the difference in S/N. 
Reassuringly, our BHB stars tend to have lower metallicities. 

8.2. Proper motion 

The proper motions of objects in SDSS are found automatically by matching the stellar po- 
sitions to those in the astrometric catalogue USNO-A V2.0 (Monet et al. 1998) derived from the 
digitized Palomar Observatory Sky Survey (POSS) plates, for objects bright enough to have a POSS 
counterpart. The normalized distribution of the proper motion of spectroscopically confirmed QSOs 
is shown in figure 12 as a thick line. These alleged proper motions of QSOs are obviously due to 
measurement errors, both from POSS and SDSS (primarily from the former). If both measurements 
follow bivariate gaussian errors with radially symmetric variances a\ and a 2 , then the distribution 
of the separation of the measurements (i.e., proper motion) also follows a bivariate gaussian with 
variance a 2 = a\ + cr|- Thus, the probability density distribution of observing a star with true 
proper motion /j,q at an observed proper motion \i is 

, , v 1 r 2 * ( +j4 - 2 wo cos9 \ n ( n 2 + nl \ /mn m 

^ W = 2^{ 6X H 2^ )» de = ^ e * P {-^^) 10 (4) 

where Iq(x) is the modified Bessel function of the first kind. Since /jlq — > for the QSO distribution, 
the second moment of the distribution of fi should equal 2a 2 for gaussian statistics. This determi- 
nation yields cqso = 4.77 mas yr -1 , and the resulting theoretical error distribution function from 
equation 4 is shown in figure 12 as a thin line. 

We have estimated the statistical intrinsic proper motions of our sample using the derived 
photometric distances. The transverse velocity of each star was drawn 1000 times at random from 
an isotropic velocity distribution with dispersion 100 km s" 1 in one dimension and corrected for 
an assumed solar velocity v = (10, 225, 7) kms -1 (Dehnen & Binney 1998). These choices were 
based on the kinematic results presented in Paper II. The distribution of intrinsic proper motions so 
derived, N(/j,q), is shown in figure 12 as the dotted line. The theoretical proper motion distribution 
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for BHB stars with measurement errors should follow the form 

POO 

<^bhb(^) oc / Ar(//o)<?Vo(>) (5) 

J 

and is drawn in figure 12 as the dashed line. The distribution of measured proper motions of 
our BHB sample is also shown. This distribution is clearly not very much broader than that 
of the QSOs; the second moment yields orhb = 4.84 mas yr -1 , which is actually too small to 
account for the combination of measurement error (ctqso = 4.77 mas yr _1 ) and intrinsic proper 
motion (<TBHBint — 2.32 mas yr" 1 )! In contrast, the proper motion distribution of a sample of faint 
(here r > 16.7) F stars from SDSS is also shown, with a much broader distribution. This finding 
demonstrates that contamination in the BHB sample from F stars must be very small. We estimate 
that the F star contamination is less than 5% among the 837 stars in our sample that were matched 
to USNO-A. Without better understanding of systematic errors, we cannot claim to have detected 
proper motion of the BHB sample. However, the USNO-B catalog, which has appreciably better 
data quality than USNO-A, will soon be provided as the first epoch for proper motions, so we defer 
further analysis to later work. 



8.3. A map of the Galactic halo and the Sagittarius dwarf tidal stream 

Figure 13 shows two views of a three-dimensional representation of the locations of the BHB 
stars. Each star is colored according to its line-of-sight radial velocity with the solar velocity 
of (10,225,7)kms- 1 (Dehnen & Binney 1998) added. These two views compose a stereograph 
which can be viewed by focusing the eyes on a point midway to the paper. The three-dimensional 
positions of the stars are derived from their angular positions and photometric distances (§7), 
assuming Rq = 8kpc. 

One could search for substructure in the Galactic halo by a careful examination of figure 13, 
and indeed a clump appears in the upper left part of the plot (at (£,b) ~ (350,50)). This clump 
belongs to the leading arm of the Sagittarius dwarf tidal stream (Ivezic et al. 2000; Yanny & 
Newberg et al. 2000; Ibata et al. 2001a, b). It is interesting to investigate this structure in velocity 
phase space, as shown in figure 14. The Sagittarius stream approximately follows the great circle 
on the celestial sphere with pole (a, 5) = (308,58) or (£,b) = (94,11), as determined from the 
distribution of SDSS candidate RR Lyrae stars. Practically the same position is implied by the 
distribution of 2MASS M giants (Majewski et al. 2003). This great circle intersects the footprint 
of the SDSS in three places. These three intersections with BHB stars (two in the north and one 
in the south, with a strip width of 10 degrees) are represented in the three panels in figure 14. One 
can see a clump in the phase space of the northern subsample representing the Sagittarius stream 
(left plot); indeed it is the same clump in the upper left of figure 13. The coherent structure in 
the middle panel (at D ~ 30kpc) is also associated with the Sgr dwarf tidal stream (trailing arm). 
This is the same clump that was discovered by Yanny & Newberg et al. (2000) using A-colored 
stars, and is also detected in the velocity-distance distribution of SDSS candidate RR Lyrae stars 
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(Ivezic et al. 2003). The fact that we see such structure in the BHB sample strongly indicates that 
our BHB absolute magnitudes and selection criteria are reliable. 

9. Summary 

Using broadband colors and spectroscopic surface-gravity indicators, we have selected 1170 
blue horizontal branch stars from the SDSS. For the 733 bright stars (g < 18) the contamination is 
probably < 10%, and for the 437 faint BHB stars (g > 18) the contamination is probably ~ 25%. 
The absolute magnitude of BHB stars was determined from globular cluster color-magnitude di- 
agrams, and an absolute-magnitude-color relation established. This was used to estimate photo- 
metric distances, with errors ~ 0.2 mag. The stars also have measured radial velocities with errors 
~ 26kms _1 . Proper motions determined by matching against the digitally scanned Palomar Sky 
Survey (USNO-A, Monet et al. 1998) also indicate that contamination from the disk is small. We 
emphasize, however, that no kinematic information was used to select the sample. The data can 
be downloaded as Table 3 in the electronic version of this paper. 

We thank A. Gould, B. Paczyhski, R. Lupton, and J. Hennawi for helpful discussions, and 
M. Strauss for a close reading of a draft of this paper. GRK is grateful for generous research support 
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the Participating Institutions, the National Aeronautics and Space Administration, the National 
Science Foundation, the U.S. Department of Energy, the Japanese Monbukagakusho, and the Max 
Planck Society. The SDSS Web site is http://www.sdss.org/. The SDSS is managed by the 
Astrophysical Research Consortium (ARC) for the Participating Institutions. The Participating 
Institutions are The University of Chicago, Fermilab, the Institute for Advanced Study, the Japan 
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Planck-Institute for Astronomy (MPIA), the Max-Planck-Institute for Astrophysics (MPA), New 
Mexico State University, the University of Pittsburgh, Princeton University, the United States 
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A. Remark on contamination and incompleteness 

There is no independently confirmed subsample of BHB stars, but we can estimate the accuracy 
of our spectroscopic selection methods by intercomparing their results. 

After applying our color cut, we are left with N + C = 2338 candidates bright enough to be 
subjected to our spectroscopic selection criteria. Of these, N are true BHB stars, and the other 
C are potential contaminants. Let rji < 1 be the probability that method i fails to identify a true 
BHB as such, and k- be the probability that it mistakenly accepts an individual member of the 
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non-BHB group. Then the expected number of stars accepted by method i is 

Ni = (1 - V i)N + n'fi 
Similarly, the numbers of stars selected jointly by methods i & j, or by i, j, & k , are 



(Al) 



A;n;nfc = (1-^(1-^(1-^ +4j4^, 

respectively, assuming that these methods are truly independent. It is convenient to define Ki = 
K^C/Ni so that the number of contaminants among the Ni stars accepted by method i is KiN. 
Then 
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As described in §4, we have four methods. Therefore we have six double intersections (A3) and 
four triple intersections (A4) with which to solve for the six unknown parameters (k±, . . . K4, N, C) 
that they contain. So the problem is overdetermined, allowing us perhaps to test the assumption 
of independence. Perfect consistency is not to be expected since even for fixed values of the above 



parameters, the numbers {N, iVjnj, ^Vjninfe 



} are subject to statistical fluctuations. 



If the contamination fractions Ki are reasonably small — it is sufficient that k/(1 
y/C/N — then the terms involving C in eqs. (A3)-(A4) are unimportant, so that 



1 



Ki 



N, 



injnk 



N: 



(A5) 



This approach yields the results shown in Table 5, and also an estimate for the true number of 
BHBs tested (i.e, not including those too faint for spectroscopic selection): N ~ 772. Note that 
there are three possible choices of (j,k) with which to estimate Ki from eq. (A5), all of which are 
shown in the Table to give an idea of their consistency. It can be seen that the contamination and 
incompleteness fractions for each are indeed small, < 10%, if these methods are indeed independent 
as we have assumed. In support of the latter, the differences among the estimates in columns 2-4 
are ~ ±0.01 as would be expected by chance since, for example, binomial statistics predict that 
the number of false negatives for method i should fluctuate by ±\/ ?7i(l — rii)N, or about 7 out of 
772 if r]i = 0.07. 
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Fig. 1. — SDSS color-color diagram showing all spectroscopically targetted objects brighter than 
g = 18 (for clarity) which were subsequently confirmed as stars. The main sequence runs from the 
center of this diagram towards the upper right, and the large Balmer jump of A-colored stars places 
them in the offshoot, where our "color cut" selection box is drawn. Note that the relative densities 
of stars in this plot are affected by selection bias; for example, more stars appear here which lie in 
the more "interesting" regions of color-color space for QSOs. 
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Fig. 2. — Spectrum of a typical (high-S/N) BHB star (upper panel) and the H7-EW region of the 
same star with the continuum divided out (lower panel). The parameters (/ m , -D0.2) are labelled 
for both lines. 
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Fig. 3. — The extracted and normalized spectrum of a BHB star (solid) and a BS star (dotted) 
clearly showing the BS star's wider Balmer lines at 80% of the continuum. 
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Fig. 4. — The parameters / m and Z?o.2 as determined by the H5 Do. 2 method for stars that pass 
the color cut with g < 18. The trail of stars with / m > .4 are too cool to be BHB stars, and the 
concentration of stars with .Do. 2 > 28 is due to blue stragglers with higher surface gravity. The 
~ 10 stars that lie well above the main locus for / m > .4 are placed there by poor parameter 
determinations due to missing spectroscopic data at the location of the H<5 line. 
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Fig. 5. — Parameter determinations of the H<5 23>o.2 method (left panels) and of the H7 scalewidth- 
shape method (right panels) for five g magnitude bins, labelled. Errors evidently become too large 
to safely discriminate BHB stars from other types of stars at about g ~ 18. 
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Fig. 6. — Comparison between the Z?o.2 method (left panels) and the scalewidth-shape method 
(right panels) for the H7 line (top panels) and the H<5 line (bottom panels). Each star is represented 
by the same color in all four plots, so the compatibility among all four methods is seen. Here stars 
are colored according to their classification by the -Do. 2 method for the H5 line: BHB stars are blue, 
BS stars are green, and the cooler stars are red. Only stars with g < 18 are plotted here. 
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Fig. 7. — Parameter determinations from the H5 Z?o.2 method for stars with g < 18. Stars are 
colored according to their location on the b — c diagram for the H7 scalewidth-shape method: stars 
for which b < 12 — 30(c — 1.18) 2 are colored blue; all others are colored red. (The temperature 
criterion is not used.) Thus stars which have low values of b from the scalewidth-shape method 
also have low values of Dq.2i and are likely horizontal branch stars. 
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Fig. 8. — The same color-color diagram presented in figure 1 but zoomed into the color cut box 
and with stars colored according to their determined classifications by the combination cut. Only 
stars with g < 18 are shown. BHB stars are blue, BS stars are green, and all others (F stars, 
cooler horizontal branch stars) are red. Note that the three samples do comprise different regions 
of color space, suggesting that purely photometric separation criteria are possible, but that the 
overlap between the regions is quite significant. The sequence of theoretical horizontal branch 
models discussed in §7 is plotted as the thick curved line. The five black circles on the thick line 
denote, starting at the top and progressing down and to the left, stars with absolute magnitudes 
g = (0.6,0.55,0.6,0.7,0.8). 
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Fig. 9. — Discrepancies between the radial velocity determinations Av r from duplicate observations 
of the same stars vs. g magnitude. The velocity error increases for faint stars, as expected. Stars 
are colored blue if both observations indicate a BHB star (using the combination cut for g < 18 
and the stringent color cut for g > 18), red if both indicate a non-BHB star, and green if the two 
observations are in conflict. The reason why there are some conflicting classifications for g > 18 
is because the two spectroscopic observations sometimes actually reference different photometric 
observations when available. 
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Fig. 10. — Color-magnitude diagrams for the Galactic globular clusters NGC 2419 and Pal 5. The 
distance modulus has been subtracted from the dereddened g magnitudes to arrive at directly 
observed absolute magnitudes. Lines at g — r = (—0.4, 0.0) denote the color cut. The theoretical 
horizontal branch sequence is overplotted. 
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Fig. 11. — Distribution of the / m parameter for the Call K line for BHB stars and non-BHB stars. 
Although BHB stars clearly have larger values of / m and thus lower metallicity implying that they 
are halo objects, they are not so cleanly separated from the non-BHB stars. This is probably 
because SDSS spectroscopy does not have high enough resolution and S/N to probe the Call K 
line as accurately as it can the Balmer H7 and H<5 lines. 
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Fig. 12. — Distributions of measured proper motion of SDSS objects, with integrals normalized 
to unity, of SDSS objects compared to USNO catalogs. The QSO histogram reflects measurement 
error. The smooth thin curve is the expected analytic form from equation 4 of the QSO distribution 
with variance (4.77 mas yr" 1 ) 2 . Our sample of 1170 BHB stars is plotted as jagged (not boxy) for 
clarity. The dotted line is the expected distribution of true proper motions for the BHB sample. 
The dashed curve is the expected measured distribution of BHBs given the theory of measurement 
errors described in the text. The close match between the QSOs and BHB stars qualitatively 
indicates that contamination from foreground objects is low. 
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Fig. 13. — Two three-dimensional projections of BHB stars discovered by the SDSS spectroscopic 
survey, comprising a stereograph which can be viewed in 3-D by focusing the eyes on a point 
midway to the paper or with a stereoscope. Each star is colored by its line-of-sight radial velocity 
with the solar velocity subtracted out. One can search for substructure in the halo by looking for 
streamers in this figure. Although substructure is not the emphasis for this work, we note the lump 
in the upper left which is the Sagittarius dwarf tidal stream. The bracket in the middle of the 
figure represents a coordinate system centered on the Galactic center with the longest line segment 
pointing towards and terminating at the location of the sun at 8 kpc. Galactic north is up. Only 
stars within a cube of half-length 50 kpc are drawn. The axes are labelled in kpc. 
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Fig. 14. — Phase space plots for three subsamples of the BHB sample extracted from the Sagittarius 
stream. The three regions of extraction are given by the intersection of the BHB sample with the 10 
degree wide strip following the great circle described in §8.3. The left plot, from the north Galactic 
cap, clearly shows a clump in phase space at ~ 50kpc which represents the Sagittarius stream. 
The middle plot, from the south, can also be seen to show similar structure. The right plot is from 
the north but the number of stars is too small to make any conclusions about the detection of the 
Sagittarius stream in this region of space. 
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Table 1. BHB criteria for the different methods 



method 


temperature criterion 


linewidth criterion 


L>0.2 method (H7) 


/m < 0.37 


D .2 < 28 A 


D .2 method (R5) 


fm < 0-34 


D0.2 < 28 A 


scalewidth-shape method (H7) 


c > 1.00 


6< 12.0-30(c-1.18) 2 


scalewidth-shape method (H<5) 


c > 0.95 


b< 12.5-40(c-1.10) 2 



Table 2. Predicted absolute magnitudes and colors for blue horizontal branch stars in the SDSS 
system, assuming eq. (3) and [M/H] = —2.0 (1st line at each T e g) and [M/H] = —1.0 (2nd line). 
In the second column g is gravity; in all other instances it is the magnitude in the g band. 





logff 


0(10 pc) 


u-g 


g-r 


9 ~ i 


7000 


2.900 


0.7115 


1.0335 


0.0945 


0.0736 






0.6862 


1.0836 


0.0999 


0.0724 


8000 


3.132 


0.5712 


1.1792 


-0.1313 


-0.2700 






0.5494 


1.1944 


-0.1291 


-0.2724 


9000 


3.336 


0.6311 


1.0695 


-0.2300 


-0.4301 






0.6121 


1.0732 


-0.2317 


-0.4351 


10000 


3.520 


0.8000 


0.8911 


-0.2742 


-0.5023 






0.7753 


0.8853 


-0.2780 


-0.5093 


11000 


3.685 


0.9954 


0.7381 


-0.2999 


-0.5449 






0.9672 


0.7239 


-0.3034 


-0.5512 


12000 


3.836 


1.1884 


0.6139 


-0.3203 


-0.5785 






1.1582 


0.5950 


-0.3232 


-0.5836 


13000 


3.975 


1.3735 


0.5103 


-0.3382 


-0.6075 






1.3425 


0.4881 


-0.3404 


-0.6115 



Table 3. List of 1170 BHB stars selected from SDSS 



astrometry uncorrected magnitudes extinction-corrected magnitudes 



ra 


dec 


e 


6 


u 


9 


r 


i 


z 


a 


9 


r 


i 


2 


redshift 


Vlos 


D (kpc) 


D 0.2,H 


145.82320 


-1.17255 


237.13901 


36.66249 


17.89 


16.68 


16.85 


16.95 


17.04 


17.65 


16.50 


16.73 


16.86 


16.97 


0.0005883 


24.71 


15.34 


27.364 


160.36295 


-1.05517 


249.66811 


47.91172 


18.08 


16.83 


16.91 


16.94 


17.06 


17.80 


16.62 


16.76 


16.83 


16.98 


0.0005751 


33.98 


16.36 


25.961, 


161.66508 


0.06255 


249.84621 


49.62371 


18.55 


17.20 


17.36 


17.47 


17.50 


18.34 


17.05 


17.25 


17.39 


17.44 


0.0005553 


32.86 


19.98 


24.1511 


161.08659 


1.00213 


248.17850 


49.86035 


18.64 


17.37 


17.40 


17.51 


17.55 


18.41 


17.20 


17.28 


17.42 


17.48 


0.0008044 


109.62 


21.27 


23.597 


161.20377 


0.68167 


248.66331 


49.72445 


17.88 


16.60 


16.72 


16.81 


16.89 


17.65 


16.43 


16.60 


16.72 


16.82 


0.0006674 


67.72 


14.98 


26.983' 



Note. — The first four columns contain the astrometry (ra, dec, I, b) for each object. The magnitudes (u, g,r,i, z) are in the next ten columns: uncorrected and corrected for extinction. 
kms" 1 , assuming vq = (10, 5, 7) kms -1 . The next column is the distance in kpc as determined by the method of §7. The last nine columns are the linewidth parameters from the Balmer an 
printed edition contains only a sample. 

I 

Co 
Co 

I 
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Table 4. Numbers of stars jointly selected 







SW-S (K'y) 


SW-S (R5) 

K-f ¥ T k_/ \ 111/ 1 


Ni Nint ■ ■ ■ 














2338 











1 


741 








1 





746 








1 


1 


683 





1 








779 





1 





1 


661 





1 


1 





658 





1 


1 


1 


620 


1 











777 


1 








1 


670 


1 





1 





692 


1 





1 


1 


642 


1 


1 








673 


1 


1 





1 


620 


1 


1 


1 





630 


1 


1 


1 


1 


595 



Note. - - The first four columns are binary bits indicating which 
selection methods were applied. The top line is the total number of 
candidate stars {g < 18); the bottom line is the number selected by 
all four methods. 



Table 5. Accuracy of four spectroscopic selection methods. 



method i 3 estimates of Ki (m) (rji) 



A).2(H7) 0.0706 0.0640 0.0486 0.0611 0.0552 

A).2(H<5) 0.0995 0.0789 0.0977 0.0920 0.0840 

SW-S(H 7 ) 0.0283 0.0131 0.0312 0.0242 0.0573 

SW-S(H5) 0.0173 0.0381 0.0360 0.0305 0.0696 



Note. - Estimates in columns 2-4 are based on the 
three triples involving i drawn from methods (1,2,3,4) 
[eq. (A5)], column 5 is their average, and column 6 is 
computed from column 5 and eq. (A2). 



